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Abstract We resume the discussion about irradiation-driven mass transfer cycles in semi-detached compact binary systems. The 
analytical model that describes the onset of these cycles, which occur on a thermal timescale of the donor star, is reexamined. We 
take into account a contribution of the thermal relaxation which is not related to the irradiation of the donor star and which was 
neglected in previous studies. Cataclysmic variables (CVs) containing extended giant donors are more stable than previously 
thought. CVs close to the upper edge of the period gap can undergo cycles for low angular momentum loss rates, as they have 
been suggested by recent magnetic braking prescriptions, while they are stable for high braking rates. 

A model for the irradiation geometry that takes into account surface elements near the terminator of the donor star indicates 
that possibly also low-mass X-ray binaries (LMXBs) can undergo mass transfer cycles. Regarding the braking rate, which is 
necessary to drive cycles, basically the same restrictions apply for short period LMXBs as they do for short period CVs. We 
confirm that LMXBs containing giants can undergo cyles. In terms of an irradiation efficiency parameter a CVs are susceptible 
to the irradiation instability for a > 0.1 while LMXBs are susceptible for a < 0.1. 

The predictions of the analytical model are checked by the first long-term evolutionary computations of systems undergoing 
mass transfer cycles with full ID stellar models. For unevolved main sequence (MS) and giant donors the analytic model 
provides reasonable values for the boundaries of the stable and unstable regions while CVs containing highly evolved MS 
donors are more stable at high braking rates than expected. 

Taking into account irradiation the minimum period of CVs is increased by up to 1-2 minutes, depending on a. 
Key words, binaries: close - novae, cataclysmic variables - stars: evolution - stars: mass-loss - X-rays: binaries 



^ 1. Introduction periodic symmetric irradiation and concluded that a more re- 

• • alistic model is required. 

_ > Observations of close binary systems have shown that a star ^ ^^^-^^ p^p^^^ iRitter et alJ (Hooj fl99^ : iKingetalJ 

^ ^ which is illuminated by its companion can exhibit a hotter, Olu- ^g^^^^ developed a feedback model for asymmetrical ir- 

■ minated, and a cooler, uniUuminated side (see, e.g. , Eddingtonl radiation that predicts under which conditions the secular mass 

.C^ ; 119261 I Ritter etalj | 2000t and references therein). This is the ^^^^^^^^ ^^^^ ^^^^^^^^ ^^^^^^^^j^ ^^^^^^ undergoes mass 

well-known reflection effect (for a review see Vaz 1985). If the ^^^^^^^^ ^^^j^^ ^j^^^j^ ^^^^^ ^ ^^^^^^^ timescale of the con- 

lUuminated star has a convective envelope, then irradiation in- ^^^^^^^ envelope of the donor star 

hibits the energy transport through the illuminated outer layers ,,, . , , . , • ,• • j- „ , 

„ ,,„n,T . .. .. We reexamine the above-mentioned irradiation feedback 
( Vaz & Nordlund 1985). In the case of a semi-detached com- 
pact binary system the reflection effect can be caused by ac- 
cretion luminosity that is released nearby the compact star In 
this case the energy transport, i.e., the intrinsic flux through 
the illuminated outer layers of the donor star is coupled to the 
mass transfer rate. This leads to a feedback, the consequences 

of which are the subject of this paper. Unfortunately, fliere have been basically no numerical 

computations to test the validity of this model and espe- 

The influ ence of irradiat iononto mass transfer was first ^.^^^^ boundaries of the stable and unstable regions. Only 

discussed by Podsiadlowskil ]1993 for the simplified_case computations have been carried out by using bipoly- 

ofs^mmetrical iiTadiation of flie donor Later, Hameury et alj j^^p^ ^j^jj^^ ^^^^^j^ ^^^^ ^^^^^ ^j^^ ^^-^^ 

119931) tried to simulate asymmetrical irradiation by time- fled models for the ti;eatment of irradiation (Ritter et al. 19951 

iMcCormick & FrankI 19981 IRitter et al1l200ol) . There has been 

Send offprint requests to: Hans Ritter, only one short-term evolution using full ID stellar models and 

e-mail: hsr@mpa-garching.mpg.de a realistic treatment of irradiation which was published by 



model and take into account a contribution of the thermal relax- 
ation which is not related to the irradiation of the donor star. It 
turns out that this term, which was neglected in previous studies 
of the above-mentioned analytical model, can become impor- 
tant, especially for giant donors. 
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iHameurv & Rittei] ( Il997l) . It has served more as a "proof of 
concept" for their model of the reduction of the intrinsic flux 
by illumination than as a numerical confirmation of the irradi- 
ation feedback model. 

The main purpose of this paper is to confirm the predictions 
of the irradiation feedback model by numerical long-term evo- 
lutionary computations that are based on full ID stellar models 
for the donor star, a realisti c model for the irradiatio n geom- 
etry, and tabulated results of IHameurv & RitteJ (Il997t) for the 
reduction of the intrinsic flux by illumination. 

In Sect. 12 we mention the input physics that is necessary 
for the feedback model and for the numerical computations. 
This includes the mass loss rate prescription, the reflection ef- 
fect, the irradiation geometry, and the thermal relaxation of the 
donor star Afterwards, in Sect.|3lwe develop two ordinary dif- 
ferential equations which describe the system adequately for 
our purpose, and in Sect. |4] we derive the conditions which de- 
scribe the onset of instability by means of linear stability anal- 
ysis. 

We give a short overview in Sect. |5] about input physics, 
which we use in our code, and subsequently, we present our 
numerical results for CVs and LMXBs containing evolved and 
unevolved donor stars in Sect. |S] 

2. Input physics 

We consider a semi-detached compact binary system. The com- 
pact primary star, a white dwarf, a neutron star, or a black hole, 
is denoted by subscript 1, the secondary star, a main sequence 
star with a convective envelope or a giant, by subscript 2. Both 
stars move on circular orbits with an orbital distance A. We as- 
sume that the secondary corotates with the orbital motion. This 
is justified by the short timescales of orbital circulatization and 
synchronization for such systems (e.g.. lzahrul977.) . The mass 
ratio of both stars is defined as 



M2 

Ml ' 



(1) 



The secondary almost fills its critical Roche lobe and therefore 
loses mass through the inner Lagrangian point Li to its compact 
companion. 

2.1. The mass transfer rate 

Since the donor star does not have a sharp rim, the mass loss 
rate must be a continuous function of the difference 



AR^Ri- Rr,2 



(2) 



between the radius R2 and the Roche radius 2 of the star. 
For computing the mass loss rate we use 



IAR\ 

M2 = -Mo expl — I, 



(3) 



where Mq > is a weakly varying function of the system pa- 
rameters Ml, M2, A, and the photospheric values of the donor, 
and Hp denotes the photospheric pressure scale height (for de- 
tails see iRitter..l988.) . Strictly speaking, Eq. (|3} is applicable 



only to a donor star which slightly underfills its Roche lobe. 
This is the case for most of the systems we are interested in. In 
any case, the conclusions of this paper are independent of the 
exact form of the mass transfer prescription as long as the char- 
acteristic scale length H, on which M2 increases by a factor of 
e, fulfiUs 



H 



1 dM2 



«/?2. 



(4) 



,M2 dAR 

We note that in (|3} the characteristic scale length is H = Hp 
and that ^ is of order ~ 10"^ for low-mass MS stars and up to 
~ 10"^ for giants. 

2.2. The accretion geometry and the irradiating flux 

We assume that the matter streaming through the Lj -point is 
collected in an accretion disk surrounding the primary and that 
a fraction < 77 < 1 of the transferred matter is finally accreted, 
i.e.: 



Ml - -rjM2- 



(5) 



The fraction 1 - 77 is lost from the system and carries away 
angular momentum, which is parametrized in terms of 

V , (6) 

dlnM ^ 

where / is the orbital angular momentum and M = Mi -1-M2 the 
total mass of the system. This mode of orbital angular momen- 
tum loss is usually referred to as consequential angular momen- 
tum loss (CAML, see, e.g. JKing & Kolbll995h . If the matter is 
lost with the specific orbital angular momentum of the primary, 
then V - q. 

The accretion process releases gravitational energy and 
produces accretion luminosity, mostly on or nearby the com- 
pact star so that the source of the accretion luminosity can be 
treated as a point source. Hence, the irradiating flux at a dis- 
tance d is given by 

rTMi 



And^ Ri 



-M, 



(7) 



r is the gravitational constant, and a.^cci takes into account that 
the gravitational energy is not necessarily completely released 
as accretion luminosity and that the accretion luminosity is not 
necessarily radiated isotropically with respect to the primary. 

On the long term the system can lose additional matter and 
angular momentum, e.g., by nova explosions so that the long- 
term avera ge value Ti can be less than ??, even ri < is p os- 
sible (e.g.. Hameur v et al.lll989l IPriahiik & Kov etz"l99^. In 
this case i) determines the (long-term) mass and angular mo- 
mentum balance while 77 determines the momentary accretion 
rate and therefore the accretion luminosity. As has been shown 
by Schenker et al. (1998), describing the long-term evolution 
by an average 77 is reasonable. 

2.3. The reflection effect and the intrinsic flux 

A star that is illuminated by its companion shows a reflection 
effect (for a review see ■Vaz,.1985i) . Accordingly, the star has 
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an illuminated, brighter and hotter side with an irradiation- 
dependend effective temperature Ti„ and an unilluminated, 
darker and cooler side with an effective temperature Tq. The in- 
trinsic flux Fi„i through a surface element is defined as the net 
flux, i.e., the difference between the incoming flux Fi„ and the 
outgoing flux err"* , where cr denotes the Stephan-Boltzmann 
constant. It is convenient to write the irradiating flux in units of 
the intrinsic flux Fq = ctTq on the unilluminated side: 



fin 



Fo 



(8) 



Typically, the reflection effect is discussed in terms of the 
bolometric reflection albedo 



o-T* - o-T^ 



(9) 



but for the purpose of this paper it is more convenient to use 
the intrinsic flux in units of Fq: 



f. - — 



err.-* 



— iri" r 

~ "Tf^ ~ /in 



(10) 



For stars in radiative equilibrium the irradiating flux fin- is re- 
processed and completely reemitted. They show Wboi = 1 and 
/int = 1 all over the surface. Stars with convective envelopes 
are by definition not in radiative equilibrium and for them vvboi 
and /int can be less than unity. 

Convection in the interior o f stars i s efficient, i.e., adiabatic. 
It was first stated by 'Rucinski' that the entropy, which 

is constant in the adiabatic convection zone, must determine 
the intrinsic flux throug h the outer layers . This means that, as 
was first recognized bv IVaz & N ordlund* ("1985'), the tempera- 
ture gradient must be flatter on the illuminated, hotter side than 
on the unilluminated, cooler side: Since the temperature gra- 
dient in the adiabatic convection zone is essentially tied to the 
adiabatic temperature gradient, the temperature gradient in the 
subphotospheric, superadiabatic layers must be flatter on the il- 
luminated side. Therefore, energy transport through the illumi- 
nated outer layers is inhibited. In fact, the superadiabatic zone 
acts like a valve for the energy flo w. This can be unterstoo d in 
ter ms of a sim ple one-zone model jRitter e t al.'l995'.'2000'). 

iRitter et al. (20D_Q. App. A) have pointed out that lateral 
energy transport within the superadiadiabatic convection zone 
is negligible since the radial temperature gradient in the su- 
peradiabatic layers is much larger than the lateral temperature 
gradient. Hence, irradiation can be treated as a local effect and 
for every surface element separately. While this is most likely 
a good approximation for weak irradiation, i.e., typical CVs, 
it is possible that for sufficiently strong irradiation, i.e., typ- 
ical LMXBs, non-local effects like circulations, which trans- 
port heat from the illuminated to the unilluminated side, be- 
come non-negligible. We will discuss this topic in more detail 
in Sec. 1641 

By solving the equations of stellar structure for an unirra- 
diated and an irradiated surface element simultaneously for a 
fixed entropy at the bottom of the adiabatic zone it is possi- 
ble to compute fnt (fn) numerically. This has been done by 



iHameurv & RitteJ ( Il997l) ' for grey atmospheres and perpen- 
dicular irradiation for a large grid of effective temperature To 
and surface gravity g - The results are available in tabular 
form and we use them for our numerical calculations. 

The assumption of a grey atmosphere and an irradiated 
Planck spectrum is not necessarily very accurate, especially 
for CVs and LMXBs. Various properties like a different pen- 
etration depth of optical light, UV and X-rays are neglected. A 
smaller penetration depth means that the external radiation field 
is mainly absorbed and reprocessed in the optically thin^ lay- 
ers of the atmosphere. This reduces the effect of the irradiation 
on the photospheric boundary conditions and therefore the in- 
trinsic flux is less inhibited. This "efficiency" of the irradiation 
for a non-grey atmosphere compared to a grey atmosphere will 
be taken into account formally by a free parameter ai^- If the 
external radiation field has the same spectrum as the irradiated 
star, then ffin is unity by definition. In the case of LMXBs hard 
X-rays, which can penetrate into the subphotospheric layers, 
can be the dominant component of the irradiating spectrum. 
But since even in the case of a grey atmosphere (i.e., ai„ - I) 
about exp^-|) ^ i of the irradiating flux penetrates directly 

down to the pho tosphere at an optical depth of r = | (see e.g., 
iTout et al.l l989. Eq. 29) and since even hard X-rays, except for 
the most energetic p hotons, do not p enetrate sufficiently deep 
into the outer layers ( lHameurvll99d) . we do not expect that ai„ 
can become significantly greater than unity in typical LMXBs. 
Therefore, in general < 1 . 

Models of non-grey atmospheres irradiated by Planck spec- 
tra (^rdlund & Vaz 1990; Brett & Smith 1993) a nd real s tellar 
spectra ( Barma n & Hauschildti2001u2002<:iBarmanl20o3) have 
been computed. However, those computations are much too 
time-consuming for generating tables that would be apphca- 
ble for our purpose, i.e., similar to those bvl Hameurv & RitteJ 
Il997). Nevertheless, it is possible to compare models of ir- 
radiated non-grey atmospheres with models of irradiated grey 
atmospheres to get an estimate for the irradiation efficiency ai„. 

Since fnt depends on the penetration depth, irradiation is 
more efficient for a smaller opacity and thus for a smaller 
metallicity. Fig. shows fnt as a function of metallicity for a 
specific binary system. Since the source radiates mainly in the 
optical, the results for grey atmospheres (crosses) do not dif- 
fer much from those for non-grey atmospheres of the irradiated 
star, for an irradiated Planck spectrum (triangles) as well as for 
a rea l stellar spectrum (diam onds). The data shown are taken 
from lNord lund & Vaz ( 1990, table 1) and are compared with 
results from .Hameurv & Ritter. Cl997i) (asterisk). 

The reduction of fnt by irradiation also depends on the 
model of convection, especiall y the mixing leng t h / as s hown in 
Fig. 121 The data are ta ken from lNordlund & Vaj(ll990t table 2) 
and .Hameurv & Ritter. (.1997.) . 



' Those authors denote /in by x, /i„t by G(x) and jjr^ by -g(x). 

' "Optically thin" means optically thin for the outgoing radiation 
field. The same atmospheric layers can be optically thick for the in- 
coming radiation field. 
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Figure 1. The intrinsic flux /im as function of metaflicity Z for 
a star with To = 4500K and an effective gravity \ogg = 4.5 
that is irradiated by fin- ~ 1.006 with an angle of incidence 
/3 = 47.9°, i.e., eff'ectively by /;„ ^ 0.674. The star is iiTadi- 
ated by its companion with T[ - 6 000 K a nd a relative radius 
of ^ ^ The data shown are taken from lHameurv & RitteJ 
( 1 19971) (asterisk) and lNordlund & Vaj ( ll99d table 1): grey at- 
mospheres (crosses), and non-grey atmospheres irradiated by 
a Planck spectrum (triangles) and a real stellar spectrum (dia- 
monds). 
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Figure 2. The intrinsic flux fim as function of the mixing length 
/ in units of the pressure scale height Hp. Otherwise see Fig.^ 



For a plane-parallel atmosphere the external radiation field 
as a function of optical depth t and angle of incidence /3 is given 
by 



/iir(T,/3) =/i,x(0,0) COS/3 exp 



COSyS 



(11) 



where /iir(0, 0) is the irradiating flux at r = ferett & Smith' 
1993). Fig. 13 shows /im as function of p. The data from 
iNor dlund & Vaz ( 1990, table 3) have been computed for 
/in- ~ 1.006 and different values of p while the data from 
iHame urv & Ritter ( 1997.) are formally valid only for perpen- 
dicular irradiation. The corresponding data points (asterisks) 
were computed for perpendicular irradiation with an "effec- 
tive" flux /in cosyS. Since Fig.|3]does not show any significant 



difference for large angles of incidence (cos /3 — > 0), our treat- 
ment of an effective flux fi„ cos /3 seems to be reasonable. 
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Figure 3. The intrinsic flux fiat as function of the angle of inci- 
dence y8. Otherwise see Fig.^ 



Since can fluctuate on short timescales and the rela- 
tion between /int and is non-linear, the effect of irradia- 
tion does not only depend on the average fin, but also on the 
amplitude and period of those fluctuations. This is taken into 
account formally by another free parameter ad ^ 1 which ad- 
dresses the duty cycle of the irradiation. For constant irradia- 
tion Q-jj is unity. For time-varying irradiation the efficiency de- 
pends on the (time-averaged) blocked part of the intrinsic flux 
/blocked := 1 - ./nt which determines the heating of the convec- 
tive envelope. Since /nt is a convex function for typical stellar 
envelopes as shown by Fig.0] i.e.. 



/nt (fll/iT,I + fl2/n-,2) ^ Al/nt (/n-,l) + fl2/nt (/n-,2) 



(12) 



for any oi, 02 > 0, oi -1-02 = 1, which is equivalent to /.'^'^ > 0, 
the efficiency and therefore aa decreases for non-constant irra- 
diation. Splitting a mean /„ into two phases of /in-,i and /in-,2 in- 
creases the resulting mean /„( and thus decreases /blocked. If the 
illumination oscillates between an "On" state with high fluxes 
(/nt ~ 0) and an "Off' state with /n ~ 0, then we get 



ad 



Ton + T off 



(13) 



where Ton and Toff denote the average duration of the "On" and 
the "Off' state, respectively. We conclude that, e.g., nova-like 
variables have o-d ~ 1 while dwarf novae can have 1 in 

case of long quiescent phases between outbursts. 

In the following we will use the symbol /„, respectively 
Firr for the effective irradiating flux 



fin 



arj 



TM, 



-M2 



instead of (0, where the total efficiency a is defined as 



a - Q'dQ'in-a'a, 



(14) 



(15) 



and Ml is replaced by (|5}. 
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Model JNo. 


1 


2 


3 


4 


5 


6 


X, 


0.71 


0.71 


0.71 


0.05 


0.36 





M2/M0 


0.3 


0.5 


0.8 


0.45 


0.60 


0.8 




0.284 


0.436 


0.694 


0.695 


0.755 


25.81 


logTo/K 


3.552 


3.590 


3.705 


3.648 


3.738 


3.582 


logL/L© 


-1.93 


-1.41 


-0.54 


-0.77 


-0.34 


2.10 


M,e/M2 


1.0 


0.23 


0.04 


0.17 


0.035 


0.63 




0.84 


1.00 


1.50 


2.2 


2.2 


41.1 


is 


-0.31 


-0.08 


1.12 


0.60 


2.4 


0.14 




0.70 


1.04 


0.90 


0.76 


1.45 


-0.2 



Table 1. Values of several stellar models. No. 1-3 are unevolved 
MS stars, 4 is a remnant of thermal timescale mass transfer of 
an evolved MS star with an initial mass of 1 .5 Mq, an initial 
central hydrogen abundance of 0.06, an initially 0.6 Mq white 
dwarf primary, t) - 0.25, an d subsequent mass loss driven 
by strong braking according to Verbunt & Zwaan ('l98l') using 
/vz - 1 . No. 5 is similar but with an initial mass of 3 Mq, an 
initial central hydrogen abundance of 0.41, an initially \AMq 
neutron star primary, and t) is determined by an Eddington ac- 
cretion rate of 2 • 10"^ Mq /yi, otherwise f] - I. No. 6 is a giant. 




fill 

Figure 4. fim for models no. 1-4 and 6 of tabled (solid, long- 
dashed, short-dashed, dash-dotted and dot-dash-dotted line, re- 
spectively). 

2.4. The irradiation geometry 

2.4.1 . The effectively blocked surface fraction 

Using the definition of fim in ilQi the net luminosity of the star 
is given by 

-2n 



Jo Jo 



/int(/in) crT^Rl sm§d§d>f. 



(16) 



but it is more convenient to write Lint as a function of R2 and 
Tq on the unirradiated side and an effectively blocked fraction 
of the surface s: 



Lint = 471(1 - s)RlcrT^ 
with 



(17) 



( I17> is a modified Stephan-Boltzmann equation for an asym- 
metrically irradiated star The star is described as if a fraction s 
of its surface were completely blocked by the irradiation and a 
fraction I - s were unirradiated. Eq. illi allows us to describe 
an asymmetrically irradiated star by a 1 D stellar model where 
the equations of stellar structure are solved "on the unillumi- 
nated side". 

Due to the axial symmetry of the external radiation field 
(I18> can be simplified: 



/mt(/i„)] sm§d§. 



(19) 



The maximum fraction of the surface that can be blocked by 
irradiation is determined by the system geometry. As can be 
seen from Fig.|5] only surface elements with < )? < J^max are 
illuminated while surface elements with i?,nax < ^ < n can not 
see the source. »?max is defined by 

R2 
A ■ 

Since /int = 1 on unirradiated surface elements and fint > on 
irradiated surface elements, (I19> yields the maximum value of 
s: 



cos §„ 



(20) 



s < 



1 r'™ 

2 Jo 



sin I? d)? = - ( 1 - cos l^max 



Obviously, we have Sniax - 2 
ters it is of order of 0.3 - 0.4. 



) -■ Smax- (21) 

< |, and for typical system parame- 




Figure 5. Irradiation geometry. A surface element E is irradi- 
ated by the source with an angle of incidence /3 to the surface 
normal. 



2.4.2. The constant flux model 

The simplest conceivable irradiation model is the constant flux 
model (Ritter et al. 1995, 2000) which assumes that /i,T has the 
constant average value 



1 aj] FMi . 

/.rr--</in->:--g^ 4— M2 



(22) 



- ^ £ [1 - An (firr)] Sin^d^d^. (18) , . [1 - f^, «/„»] 



over the irradiated surface where (/in) denotes the irradiating 
flux at a distance A according to (I14> . Then ( I19t yields: 



(23) 
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2.4.3. The point source model 

A m ore realistic model is the point source model jKing et alJ 
Il99 6: Ritter et al. 2000) where the irradiation source is treated 
as a point source as shown in Fig. |5] In this case the effective 
flux for a surface element E is given by 



where cos fi can be expressed by 
cos & - cos 



(24) 



cos/? = 



=: h{ff). (25) 



(1 - 2 cos )5- cos )?,„ax +COS^)?max)" 

2.5. Thermal relaxation 

At turn-on of mass transfer and thus of irradiation a certain 
effective fraction s of the surface of the donor star is immedi- 
ately blocked by irradiation. As a consequence the star is out 
of thermal equilibrium because the (unchanged) nuclear energy 
production Lnuc = is higher than the (reduced) intrinsic lu- 
minosity Lint = (1 - s)Lq. The reaction of a star to sponta- 
neous blocking of its intrinsic luminosity has been discussed by 
ISpruitCl982) in connection with star spots: On short timescales 
the changes of the photospheric values on the unilluminated 
surface are negligible; on a thermal timescale of the convective 
envelope Tee the star relaxes thermally according to the modi- 
fied outer bo u ndary condi t ions b y adjusting R2, Tq, and Lnuc- 
iRitteJ (I1994I) : iRitter etal.1 (I2OO0I) give an analytical estimate 
for the effect of irradiation on the thermal equilibrium values 
of low-mass MS stars. The effective temperature remains basi- 
cally unchanged and the radius increases by about (1 - s)"° 
i.e., < 5% for typical values of s. 

For our numerical computations we compute R2 by solving 
the equations of stellar structure for the donor star. This will be 
described in more detail in Sect.|5l 

2.6. Mass transfer cycles 

Fig. |6l shows schematically how the feedback between mass 
transfer rate, irradiation, and thermal relaxation works: As soon 
as mass transfer starts, the accretion rate Mi and also the ac- 
cretion luminosity Laeer increase. This enhances the irradiating 
flux Fin- which reduces the intrinsic luminosity Lint. Thus, the 
thermal expansion rate ^2 of the star is increased which in turn 
leads to an increasing mass transfer rate -M2. 



increases 

-M2 > Ml 



L;ic 



increases decreases 

R? < Lin, < ^ir 



Figure 6. Schematic diagram of the irradiation feedback mech- 
anism at turn on of mass transfer 

Yet, this feedback can not work forever because thermal 
relaxation saturates at last after about a Kelvin-Helmholtz 



timescale tkh- But before the star reaches a new equilibrium 
state, the expansion rate starts dropping and the mass transfer 
rate is reduced. Then the feedback works into the opposite di- 
rection: Ml and Lacci decrease and so does Fin-. Therefore, Lint 
increases and the star starts to shrink. This reduces the mass 
transfer even more. In extreme cases the mass transfer even 
stops and the system appears as a detached system until the 
next cycle starts. 

3. The analytical model 

We will now derive an analytical model that describes the onset 
of mass_transfer cycles. In principle, this has already been done 
by iKing et alJ (Il996) using a slightly different formalism. AU 
expressions derived in section |3l and also 0] will be applied to 
the analytical model only and will not be used for the numerical 
computations. 

3.1. The governing differential equations 

The time evolution of the Roche radius Rr 2 can be written as 



dln/?R,2 dlnM2 -/dln7 
■ : + 2 



df 



dt 



df 



(26) 



where subscript sys denotes the systemic angular momentum 
loss at constant total mass M (e.g., angular momentum loss by 
gravitational radiation) and 

= (1 - f])(2v + 1) - 2 (1 - ^7<?) + (1 + Tjq) 4^ (27) 



I + q 



dln<7 



the niass r adius exponent of the critical Roche radius (e.g., 
lRitteiil999 )^ Angular momentum loss by CAML is taken into 
account by ^r, and the dimensionless Roche radius 



h- — 



(28) 



is available in tabulated form ("Moc hnacki' 1 9 8 4*) . 

The time evolution of the radius Ri can be splitted into three 
terms: The adiabatic reaction at constant entropy and chemical 
composition, the thermal relaxation at constant chemical com- 
position, and the nuclear evolution: 



dln/?2 dlnM, IdlnRA /din/?, 
— ^+1 — ^1 +' 



df 



df 



df 



df 



(29) 



with the adiabatic mass radius exponent fs- 

The time evolution of the equilibrium radius /?e,2 of the 
unilluminated star is given by 



dln/?e2 dlnM2 /dln/?e2 



df 



df 



df 



(30) 



wher e denotes the thermal equilibrium mass radius exponent 
re.g.. Web binkll985l) . 

With the timescale of nuclear evolution 



IdXnRiX ' /dln/?e,2y 

\ I nuc \ In 



(31) 



In that paper q is defined as i.e., inverse of the definition here. 
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the timescale of systemic angular momentum loss 

-I 



Tj 



'dln7 



sys 



and the timescale of thermal relaxation 

we define the driving timescale for the mass loss by 
1 _ 1 2 

T d Tnuc Tj ' 

and the driving timescale including the thermal relaxation by 

(35) 



(32) 



(33) 



(34) 



1 _ 1 1 

t'^ Td T"th 



Using j26l l and i29\ and taking into account R2 ~ ^r,2 the 
time evolution of AR = 7?2 - /^r 2 can be written as 



-AR = (Cs - Cv)R2— + — + — . 

df M2 T-th Td 



(36) 



Using M9\ and (13 0> and assuming^ R2 ~ /?e,2 the time evo- 
lution of 



A/?e - Rl — Re,2 

can be written as 
d M2 Rl 

df Ml Tth 



(37) 



(38) 



J36> and ( I38> are the governing d ifferential e quations and are 
analogous toEqs. [11] and [12] of lKing et ai]il996) . 



3.2. The stationary state 

One can easily see from (|3} that the mass transfer rate is sta- 
tionary if and only if AR is constant. Setting the left-hand side 
of i36\ to zero yields the stationary mass transfer rate 

M2 / 1 1 \ 1 M2 



Mi = — 



(39) 



- \ Tth Td / (s- T^ 

Additionally, the thermal relaxation 

K:=— (40) 

Tth 

is stationary if AR^ is also constant. Setting the left-hand side of 
J38> to zero and inserting ( I39> yields for the stationary thermal 
relaxation K: 



K 



1 



(41) 



Rl Tth - Td 

From ( I39> and i41\ it can be seen immediatly that the station- 
ary solution is unique if there is any at all because M2 and K 
depend on the system parameters ^s, ^e, (r, and Td only. By 
inserting ( I41> into ( I39> the timescale of mass transfer can be 
expressed by the driving timescale: 
Ml 

Tm:=-^ =(^e-^R)Td. (42) 

Ml 

This assumption is not necessarily fulfilled, especially for giants, 
but it simplifies the following equations and does not affect the results 
significantly. 



4. Linear stability analysis 

To determine the conditions for stability of the stationary 
solution we perform a linear stability analysis (see, e.g., 
iGuckenheimer & HolmesI 19831) on the vector field 



AR 

AR, 



((^s-(R)Ri^^+K+li 



{C,-C.)Ri^ + K 



(43) 



around the stationary solution ^ AR, AR^, which is a fixed point 
ofF, i.e., f (A^,A^e) = 0. 

4.1. The necessary criterion 

According to the theorem of Hartmann-Grobmann a fixed point 
X of a. vector field F is stable, if all eigenvalues of the Jacobi 
matrix DF of f at jc have a negative real part, and x is unstable 
if at least one eigenvalue has a positive real part. 

Using Q and neglecting the derivatives of Tnuc and Ri, 
which is reasonable for t„uc, tm s> tkh, and neglecting also 
the derivatives of tj we get 



DF\^ 

AR, 



Vts ^,R; Hf Ml (JAR i?ARe 



(44) 



at the fixed point. 

DF has the eigenvalues 



^1 2 = -tr ± -1/ -tr ^DF - det DF . 

• 2 V 4 



(45) 



Since DF is a real matrix, the necessary condition for 
Re/li 2 < is 



Rl — dK dK 
XiDF = - (r) —Ml + + < 0. 



(46) 



If the thermal relaxation K is neglected, Eq. J46> is equivalent 
to 

- ^R > 0, (47) 

the w ell-known criterion for dynamical stability (e.g.. lWebbinkl 
Il985l) . Since only dynamically stable systems can undergo long 
phases of mass transfer, we consider only such systems. 

The first term in ( I46> reflects the adiabatic reaction of the 
star and its Roche radius on mass loss. For dynamically stable 
systems this term is negative and therefore it contributes to the 
stabilization of the fixed point. The third term in (I46> tells how 
thermal relaxation changes with the radius of the star compared 
to its equilibrium radius. A negative value means that the larger 
(smaller) Ri, the stronger the thermal relaxation that decreases 
(increases) the radius to its (irradiation-dependend) equilibrium 
value. Therefore, 



dK 
dAR, 



<0 



(48) 



is a reasonable assumption about the thermal relaxation K. 
Thus, the necessary criterion for stability ( I46> can only be vio- 
lated if 
dK 



dAR 



>0. 



(49) 
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This means that an increasing (decreasing) AR and therefore 
an increasing (decreasing) mass transfer rate tends to enhance 
thermal relaxation which drives the expansion (contraction) of 
the star even more. Exactly this kind of feedback is caused by 
irradiation. The question is now whether this feedback can be 
strong enough to overwhelm the stabilizing terms in ( I46> . 

4.2. The sufficient criterion 

The sufficient criterion for Re A12 < is that the absolute value 
of the real part of the square root in ( 145 > is less than the absolute 
value of jtr DF, otherwise at least one eigenvalue has a posi- 
tive real part. Using (I44t it can be shown that this is equivalent 
to det DF > which is equivalent to 



- > 0. 



(50) 



This is the we ll-known criterion for thermal stability (e.g., 
IWebbinkll985l) . Only thermally stable systems can have stable 
stationary solutions since thermally unstable systems transfer 
mass on a thermal timescale. As a consequence, the thermal 
relaxation K can never reach a stationary value. Therefore, we 
consider only systems that are thermally stable. 

4.3. The onset of instability 

There are different possibilities for the dynamics of the system 
nearby the fixed point depending on tr DF and det DF if i47l 
and (15 Ot are fulfilled: 



- tiDF < -2 VdetZ>f: We have ti^DF > 4detDf so that 
the eigenvalues Ai^2 in ( I45t are real and negative. All so- 
lutions near the fixed point converge directly to the fixed 
point without any oscillation. 

2 VdetDf < trDf < 0: We have tr^DF < 4 det DF so 

that the square root in ( I45> is imaginary. The eigenvalues 
can be written as A12 = -A + ioj with A,a> > 0. All solutions 
near the fixed point spiral in. 

- ti DF = 0: The eigenvalues Ai^2 - ±i<y are imaginary. 
The system undergoes a Hopf bifurcation, when the fixed 
point, that is stable for tr DF < 0, becomes unstable, and 
a stable limit circle arises for tr DF > 0. In principle, this 
could be proved mathe matically i n ter ms of normal form 
theorv dOuckenheimer & HolmesI 19831 theorem 3.4.2 and 
eq. 3.4.11) but it requires explicit knowledge of M2 (AR) 
andK(AR,AR^). 

- < tvDF < 2 VdetDF: We have tv^DF < 4det DF so 
that the square root in ( I45> is imaginary. The eigenvalues 
can be written as Ai,2 - +A + io) with A,<jj > 0. All so- 
lutions near the fixed point spiral out and converge to the 
limit cycle sketched in Sect. 12.61 

- 2 y/det DF < tr DF: We have tr -Z)f > 4 det DF so that the 
eigenvalues Ai,2 in ( 145 > are real and positive. All solutions 
near the fixed point diverge directly without any oscillation 
and converge to the limit cycle. 



4.4. The irradiation feecibaci< 

Since K depends on AR via the intrinsic luminosity Lint, the 



derivative ^ can be written as 



dK 



dAR \ dL 



dK 



AR, 



[dAR j 



(51) 



AS, 



Lint does not depend on AR directly but indirectly via s. 
Additionally, R2 and Tq do not depend on M2. Therefore we 
can write: 



\dA/?/ 



dU 



ds dM2 



ds jj^^ j^ dM2 dA/? 



From \n\ we obtain 



diini 

d? , 
and from (|5} 



-4ncTTt,Rl 



dM2 



M2 



dAR 

With the definition of 

, ■ ds ds 
s := M2- 



'dM2 



dln(-M2) 



>0, 



(52) 



(53) 



(54) 



(55) 



which is a measure for how strongly the blocking of the intrin- 
sic luminosity changes with M2 (or /in-), Eqs. ( I51> -( I55> give 



dK 
dAR 



Hp 



dK 



dL 



int/AR, 



Using (|5|i and ( I24t in il9i leads to 

^ Jo fl/iiT 



for the point source model while ( I23> leads to 



, __1 



dA, 
'dA- 



{fin) 



(56) 



(57) 



(58) 



for the constant flux model. Figs.0and|8]show and s'^^ for 
different stellar models and also for a simple prescription for s 
that was used bv King et aLc 1995. .1997.) : 



Sk = imax tanh (k (fi^)) . 



(59) 



A: is a free parameter that defines the position of the maximum 
of 



COSh^(k{fin)) 



(60) 



As can be easily seen from ( 15 8t and also from ( I57> or i60\ . 

s' vanishes for (/„-) — > (no irradiation means no feedback) 
and for (/„-) — > 00 (strong irradiation means saturation, i.e., 
■s - Smax is constant). 
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0.20 



0.15 



0.10 



0.05 



0.00 




10^ 



Figure?. s'^_^ (Eq. ( I58» for stellar models no. 1-4 and 6 of ta- 
ble [2 (solid, long-dashed, short-dashed, dash-dotted and dot- 
dash-dotted line, respectively) and also s'j^ for k - I (Eq. i60\ . 
dotted line). 



0.20 



0.15 



0.10 



0.05 - 



0.00 




Figure 8. s'^, (Eq. (|53) for the same stellar models as in Fig.0 
For comparison of the point source with the constant flux 



model also s', for k = 1 (Eq. (I60» is shown. 



4.5. The bipolytrope model 

In terms of a bipolytrope model the thermal relaxation K can be 
expressed explicitely (Kolb & Ritter 1992, eq. 32). If the star is 
either fully convective or the nuclear luminosity of the radiative 
core is identical to the total luminosity of the core, which is the 
case in the stationary (i.e., quasi-equilibrium) state, then the 
thermal relaxation can be expressed as 



K 



Mln/^a 
df 



= -r 



TMl 



(61) 



T is given by Eq. [47] o f lKinget alJlll996ll for chemically ho- 
mogeneous stars and can be approximated by 



3Mee 



(62) 



where M ce. denotes the mass of the convective envelope 
jKing et al.. 1996; .Ritter et al..2000i) . 



The only quantity in ( I61> which depends on Lint is Lint it- 



self: 
dK 



Ro 



dLint/AR, \dLint^2 



d K 



TRl 



Inserting this into j56> gives 
dK 



dAR 



A-no-Tf.R'l 
HvYMl 



(63) 



(64) 



This can be expressed by quantities of the unilluminated star: 
In terms of the Kelvin-Helmholtz timescale 



Re,2^0 



(65) 



the corresponding timescale of the convective envelope is given 
by 



_ Tkh 3 Mce 



(66) 



(67) 



Using Hit and ( I66t Eq. i64l can be written as 

dK _ R2 s' IR, 

dAR ~ lh~,\R^ 

and is analogous to Eq. [44], respectively [59] of iKing et alJ 
(Il996l) . Alternatively, (I64> can also be expressed by quantities 
of the illuminated star: 

dK R2 s' 1 



dAR Hp I — S Tce.s 



(68) 



In this case the (irradiation-dependent) timescale of the con- 
vective envelope 



^ce,s ~ 



is defined via the Kelvin-Helmholtz timescale 



T"KH,.s 



RlL'ml 



(69) 



(70) 



of the illuminated star taking into account that it can radiate 
only through its unilluminated surface with Lint given by Ml\ . 

4.6. The homology model 

The energy generation rate from hydrogen burning can be ap- 
proximated by a simple power law of density p and tempera- 
ture T: 



pT". 



(71) 



For low-mass MS stars, whic h are considered here, the pp- 
chain gives n = 5 - 6 (see, e.g..*Kippe nhahn & Weigert.l990i) . 
By using homology relations we obtain 



R2_ 

R^,i 



-(«+3) 



(72) 



for main sequence stars. Since for giants the nuclear luminosity 
does not depend on the stellar radius, the following expressions 
also apply to giants if we set n — -3. 
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Inserting ([TtJ and (|23 into (l6ll 



gives: 



4.7. Predictions of tlie analytical model 



(n+l) 



(73) 



The equilibrium radius /?e,2 of the unilluminated star neither 
depends on AR nor on A/?e- This means 



■ (/?e,2 + A/?e) = 1- 



dRi _ d 
dAR^ ~ dAR, 

Therefore, the derivative of K with respect to AR^. yields: 



(74) 



dK 
dAR^ 



TLoR 



e,2 



4 +(«+!) (fey 



<n+2) 



(75) 



This can be written in quantities of the unilluminated star: 

\3 , / D, \-(«+2)" 



dK 
dAR. 



This result is analogous to Eq. [43] o f lKing etalJ (11996*) who 
have computed the derivative^ of ^ while we use the derivative 
of K. This leads to different values for the powers of in their 
Eq. [43] compared to our Eq. j76> : 3 and n+2 versus 4 and« + l. 
The origin of this difference is the approximate transition from 
linear to logarithmic derivatives in Eq. [2] of King et al. (1996). 

Alternatively, ( I75t can also be expressed in quantities of 
the illuminated star Using ( I69I I. ( I70> . ( I72t and Lint - L„uc in 
the (uTadiation dependend) equilibrium state we get the more 
compact expression: 



dK 



n + 5 



(77) 



Finally, we can insert (I68> and ( I77> into the necessary criterion 
for stability ( I46> and get 



dln(l - s) 
dhi(-M) 



s' Tce.s Hp 

1 < — + — (n + 5) 

l-S T . R2 



(78) 



in quantities of the illuminated star. By inserting ( I67> and ( I76> 
into ( I46> we get the equivalent result in quantities of the unil- 
luminated star: 



Trp Hp 

s' <-^ + —6, 
t' R2 



where 6 is given by 



(79) 



-(n+3) 



(80) 



If s' is less than the ratio of Tee and plus a term of the order 
of several then the stationary mass transfer rate is stable. 
Otherwise it is unstable and mass transfer cycles occur. 



Taking (|55J, eliminating A by M2 by (j^ and using the 
driving timescale including the thermal relaxation from ( I35> 
we obtain 



1 



arj 



^flTMiMi 



{fin). 



With ( I66t the ratio of Tee and can be written as 



Tee 



^s-^R^2^l4,2 



</in-). 



This result can be approximated by using \62\ . R2 
the approximation ( Pa czvhskL1971.) 



/r (?) ~ I ^ 



q + 1 



as 



(76) 1^ ^ _ _ 



3 n^^^ 
7 



ff;; M2 /(2 



(81) 

(82) 
/;r,2 and 

(83) 
(84) 



This is a function of (/!„ ) and its slope basically depends only 
on ar], the relative mass of the convective envelope jg- and the 
ratio of the radii of both stars ^ . The other factors are of order 
unity for typical system parameters. 

For the moment we neglect the last term of order of several 

in criterion ( I79> as King et al. ( 1996) did. To get mass trans- 
fer cycles the blocking of the intrinsic flux (s') as a reaction 
upon an enhanced mass transfer rate must increase faster than 
the convective envelope can relax thermally (Tee) compared to 
the driving time scale of the mass loss (r^). Therefore, systems 
are the more susceptible to the onset of mass transfer cycles, 
the smaller the timescale of the convective envelope Tee and the 
larger the driving timescale Td ~ t'^ is. 

Accordingly, main sequence stars with thin outer convec- 
tion zone s (M2 ~ 1 Mn) h a ve been sup posed to be most sus- 
ceptible dKingetalJ Il996t iRitter et al.1 12OO0I) . Also systems 
driven by the longest possible timescale, i.e., the timescale of 
angular momentum loss by gravitational radiation, have been 
supposed to be more susceptible than other systems. But there 
are basically no systems except, e.g., low-mass CVs in or be- 
low the period gap that are driven mainly by gravitational brak- 
ing. Higher mass systems are believed to be driven mainly by 
magnetic breaking. For a high braking rate, e.g., according to 
Iverbunt & Zwaan (1981), CVs with unevolved main sequence 
stars become unsusceptible to the irradiation instability except 
for the most massive systems because Td becomes too small in 
comparison to Tee (Ritter etal. 1995, 1996, 200(1). For giants 
magnetic braking is believed to be ineff'ective and Tee is small 
compared to Td ~ Tnue- Therefore, such systems and also CVs 
with highly evolved MS stars, which have a short Tee, have been 
supposed to be a priori more susceptible to the irradiation in - 
stability than CVs with unevolved MS stars ('King et alJl'9971) . 

But the last term of order of several ^ in (17^ can be ne- 
glected only if 
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As soon as Tee gets small (as for main sequence stars of ~ IMq) 
or is no longer very small (as for giants), this term must be 
taken into account. Obviously, the range of (/;„) in which 



s > — 6 

Ri 



(86) 



can be fulfilled, is an interval [(/in-.min) , (/in-,max)] because s' 
vanishes for (/in) — » and { fi„) — > oo. Thus, mass transfer 
cycles can only occur in a limited range of (/irr). This range, 
that is predicted by the analytical model, has to be checked by 
numerical computations with full stellar models. 

In this section we have derived a more concise and, in some 
sense, more general, expression ( I78> . respectively ( 17 9> for the 
onset of inst ability than the co rresponding Eq. [45], respec- 
tively [60] of'Kin g et al.l ( Il99 ^. We have also provided an ex- 
plicit formula j84> for ^ which we will use in Sect.|5]to discuss 
our numerical results. An important result is that the additional 
term of several ^ in J79> provides an upper limit for ra, i.e., 

a lower limit for -M2. Hence, for a sufficiently low secular 
mass transfer rate every system becomes stable. We will show 
in Sect. l6.1.3l that this limit can become important for CVs with 
giant donors. 

5. The binary evolutionary code 

We u s e the s te llar e volutionary code of Schlat tl et alJ 
('1997); Schlattl (1999) which goes back to the code of 
fKip penhahn et al. (1967). The c ode us es a special grid point 
al gor ithm ( Waeenhuber & Weis j fl994l). nuclear r eaction rates 
of lCau^hlan et al.j ( 1985): Adelberger et al. ( 1998), the mixing 
length theory of "Bohm- Vitense ( 1958); Cox & Giuli (.1968) . 
the OPAL opacities ( Iglesias & R ogers 19 96), and for the outer 
layers opacities by [Alexander & FergusonI ( 119941) and addi- 
tional data from P. H. Hauschildt and J. W. Ferguson (priv. 
comm.). In order to avoid numerical instabilities in calculat- 
ing mass transfer it was necessary to implement an equation of 
state which is smooth in the independent variables P, T, and 
chemical composition over the whole ra nge of application. In 
practice, we use the equation of state of ISaumon et alJ ( 1 19951) 
with data provided by I. Baraffe (p riv. co mm.), and particularly 
the equation of state of Pols et al 1 (ll995l) . Furthermore, we de- 
veloped an implicit algorithm for the trea tment of mass trans- 
fer sim ilar to the method described by Benvenuto & de Vitol 
(l2003h . For more details about our code see Buning (200 3b. For 
determining /im (fin) we use results given bv iHameurv & RitteJ 
(1 19971) and additional data from J.-M. Hameury (priv. comm.). 

6. Numerical results 

To destabilize the stationary mass transfer the ratio of Tee and 
has to be sufficiently small, at least less than the maximum 
of s' which is ~ 0. 1 . As can be seen from (I84> . especially for re- 
alistic values of a;; < 1 the ratio has to be sufficiently small 
for typical system parameters which is basically the case for 
compact binaries only. Since the system has to be dynamically 
and thermally stable and since the donor star has to maintain 
a deep convective envelope, essentially only CVs and LMXBs 
are eligible. 



Model No. 


M,/Mo 


q 


1 


V 


(r 


a 


1 


0.8 


I 


1 




-0.35 


0.3 


2 


1.2 


T2 







-1.20 


0.3 


3 


0.8 


5 
8 







-0.89 


0.3 


4 


0.8 


5 

8 





9 


-0.89 


0.1 



Table 2. System parameters of CVs containing an 0.5 Mq 
donor as shown in Fig.|5] 



6.1. Cataclysmic variables 
6.1 .1 . Unevolved MS donors 

Exemplarily for a wide range of CVs with unevolved donor 
stars Fig. |5] shows the left-hand and also the right-hand side of 
( I79> . namely s' and 



f :- 6 — -I- — 

^ Ri T 



(87) 



for a CV with an 0.8 Mq white dwarf and an unevolved 0.5 Mq 
MS secondary using system parameters listed in table|2l Low- 
mass MS stars have |f ~ 10""^ and 5 ~ n + 5 ~ 10. This 
leads to an almost constant contribution ~ lO"-' to f while ^ 
is almost linear in (/in). The most significant uncertainty in 
the values of / is due to the undetermined value of a (short- 
dashed vs. dash-dotted line in Fig.|9} while the influence of q 
(long-dashed vs. dash-dotted line) is negligible. Only 77 (dot- 
dash-dotted vs. dash-dotted fine) has a noticable influence on 
/ in this case since - is much smaller in the conservative 
case. Systems with non-conservative mass transfer (^ < 1 and 
V = g') are systematically more stable than systems with con- 
servative mass transfer (77 = 1). In the following discussion we 
will restrict ourselves to ?? - 1 and fi - Q since this case is more 
realistic than fj - \ dHameurv et aljfl989[ IPrialnik & Kovet j 
|1995). In general, the amplitude of the mass transfer cycles in- 
creases with the difference between s' and / for any given (/m). 

As can be seen from Fig. |9] CVs with unevolved donors 
are the more stable, the smaller the irradiation efficiency a is. 
Such systems can undergo mass transfer cycles if or > 0.1. For, 
e.g., a - 0.3 the "island of instabiUty", which is defined by 
s' > f, ranges from </ii-r) » 5 ■ 10"^ to </;„■> a; 2. Using (|84j is 
it possible to compute the coiTesponding ratio for Tee and t^: 



1.5- 10"^ < — <7 - 10 



-2 



(88) 



The timescale ratio must be within this range to destabilize the 
stationary mass transfer. With 



(s - 4R 



(89) 



which can be obtained from (I35> and (I41> . and Tee ~ 4.5 • 10^ yr 

for an unevolved 0.5 Mq MS star, we obtain a condition for the 
driving timescale t^'- 



3 ■ 10** yr < Td < 10"yr- 



(90) 
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10"^ 10"2 10"' lo" 10 



<fnr> 

Figure 9. s' for a CV with an unevolved 0.5 Mq MS star (no. 2 
of table[3 and an 0.8 Mq white dwarf for the point source and 
the constant flux model (solid and dotted line), s' is compared 
to / for models 1-4 listed in table |2l (dot-dash-dotted, long- 
dashed, dash-dotted and short-dashed line, respectively). The 
mass transfer is unstable and undergoes cycles if s' > f which 
is only the case in a certain interval [{/in-,min) , (/iiT.max)] whose 
limits depend on a. Otherwise, the system is stable. 




12 3 4 



P 

Figure 10. Mass transfer rate -M2 [A^o/yr] from a slightly 
evolved 0.5 Mq MS star with an age of ~ 10'° yr and » 0.62 
onto an 0.8 Mq white dwarf for the point source model as func- 
tion of orbital period P [hr]. System parameters: i] = 0, v - q, 
a = 0.3, and J = 57giav (solid line). The dotted line shows the 
mass transfer rate for the same system without irradiation feed- 
back. Both computations have been stopped beyond the mini- 
mum period at M2 = 0.04 Mq. 



If the evolution of the system is driven by gravitational braking 
only, i.e., T"d = -y ~ 4 ■ 10^ yr, then the system must undergo 
mass transfer cycles. Gravitational braking provides a physical 
upper limit for the driving timescale. Since the lower limit in 
Il90> is given by about -j^ of the gravitational braking timescale, 
we can expect the occurrence of the irradiation instability be- 
tween J = 7grav and J 



lO/grav which is roughly confirmed by 



our numerical models. 

Fig. [TO] shows a CV evolution for initial parameters of 
model no. 3 of table |2l using the point source model and an 
ad hoc braking rate of J - 57grav The donor, a slightly evolved 
MS star, fills its Roche lobe at an orbital period of about 4 hr 
and starts to transfer mass undergoing cycles. At about 3.2 hr 
the system leaves the "island of instability" and becomes sta- 
ble. The small spike in the mass transfer rate at about 2.8 hr is 
caused by the star becoming fully convective. For gravitational 
braking only as shown in Fig. ^2 the phase of mass transfer 
cycles does not end until the system is in the middle of the 
period gap. The transition, when the star becomes fully con- 
vective, coincides with a drop in the peak mass transfer rate in 
this case. 

The secular mass transfer rate is lower when taking irradi- 
ation into account (dotted versus solid line in Figs.llOland ll 1> 
since the irradiated star has a systematically larger radius. As 
a consequence, the minimum period increases, in the case of 
pure gravitational braking by up to 2 minutes, depending on 
a. This is of the same order as the correction which is caused 
by geome trical distortion of the mass losing star by the Roche 
potential dRenvoize et aljE002h . and this is still not sufficient 
to e xplain the observed orbital minimum period of ~ 78 min 
jRitter & Kolb.2003t) if only gravitational braking is taken into 
account. 

Fig.ll2lshows the phase of mass transfer cycles as function 
of time for the same system as in Fig. [TO] Since the eigenvalues 




Figure 11. Same as Fig.fTOjbut for gravitational braking only. 

at the fixed point are of the form A + ico, the frequency ai of the 
small amplitude oscillations at the end of the unstable phase 
can be computed from linear stability analysis using ( I44> and 
(|77j: 



For the system considered we have ^ 

~ Td. The period of the low amplitude mass transfer cycles 
is therefore given by 

27r 



10- 



(91) 



10 and 



r = — X 0.02 YT^^. 



(92) 



For higher braking rates, e.g., according to 
IVerbunt & ZwaanI ( Il98ll) . mass transfer is stable. Only 
CVs with more massive donor stars can undergo mass transfer 
cycles for such a high braking rate. Fig.[n]shows the evolution 
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Figure 12. Mass transfer rate -M2 [MQ/yr] as function of time 
t [yr] for the same system as in Fig.^| The computation stops 
atM2 ~ 0.34 Mq. 

of a CV with an initially 1 Mq MS donor The stationary mass 
transfer rate becomes stable when the mass of the convective 
envelope Mcs and therefore r^e becomes large enough (see also 
[Ritter et al..20na). 




P 



Figure 13. Mass transfer rate -M2 [MQ/yr] from a slightly 
evolved 1 Mq MS star (X^ ~ 0.62) onto a 0.8 Mq white dwarf 
for the point source model as function of orbital period P [hr] . 
System parameters: ij - Q, v - q, a - 0.3, and strong braking 
according to Verbunt & Zwaan (1981.). The dotted Une shows 
the evolution without irradiation feedback. 

The dependence of the temporal evolution of mass transfer 
cycles on the efficiency parameter a is shown in fig. The 
duration of the high state decreases with decreasing a until the 
mass transfer rate becomes sinusoidal with a frequency oj as 
given in (I91> . For even weaker feedback the oscillations are 
damped or even vanish completely. 



Figure 14. Mass transfer rate -M2 [Mq / yr] from an unevolved 
0.5 Mq MS star onto an 0.8 Mq white dwarf for conservative 
mass transfer (jj - 1), J - Jgav, and the constant flux model 
for a = 0.02, 0.04, 0.1, 0.2° 0.4, and 1.0 (dotted, dot-dash- 
dotted, dash-dotted, short-dashed, long-dashed, and solid line, 
respectively) as function of time t [yr] . 



counting for evolved low-mass donors in CVs is prior thermal 
timescale mass transfer from an evolved, initially more mas- 
sive star. In this case the mass losing star is more massive to 
such an extent that the system is initially thermally unstable 
(Cj-i^fi < 0) and th e mass transfer occurs on a thermal timescale 
jKolb et al.l200(ilSchenker et all2002l). 

Although the bipolytrope model is formally valid only for 
chemically homogeneous stars we tentatively apply it also to 
chemically evolved donor stars. As an example fig.^lshows s' 
and / for an 0.45 Mq remnant of thermal timescale mass trans- 
fer with a central hydrogen abundance of - 0.05. It is the 
core of an initially more massive MS star of 1 .6 Mq. According 
to (I66> Tee is less than -j^j of Tee of an unevolved MS star of the 
same mass, and its radius is larger by ~ 70%. Due to its deeper 
superadiabatic convection zone the maximum of s' is shifted 
to higher values of (/„■) and due to its larger radius / is also 
shifted to larger values of (/iir), but in total the difference to 
the case with an unevolved donor discussed before is compara- 
tively small. As can be seen from fig. El this particular system 
should be stable for a = 0.1 but can undergo cycles for a = 0.3. 
When going to smaller donor masses or lower braking rates 
during the initial evolution also remnants of thermal timescale 
mass transfer can become unstable for a > 0.1 as CVs with 
unevolved donors do. It seems that CVs with an evolved donor 
are not more susceptible to mass transfer cycles than CVs with 
unevolved donors (if at all, it is the other way around). Rather 
they are susceptible at a shorter driving timescale Td. 

As for unevolved CVs it is also possible to determine the 
range of ra where mass transfer cycles should occur, e.g., if 
Q' = 0.3: 



6.1 .2. Highly evolved MS donors 

The timescale of central hydrogen burning for low-mass MS 
stars is much longer than a Hubble time. One way of ac- 



3.5- 10^yr<Td <3.5- 10''yr. (93) 

According to this result such a system should undergo mass 
transfer cycles even for rather high braking rates. However, 



14 



Andreas Biining and Hans Ritter: Long-term evolution of compact binaries . 



- 10' 




10" 



10"' 



10"^ 



10" 




10 



20 



Figure 15. s' for a CV with an highly evolved 0.45 Mq rem- 
nant of thermal timescale mass transfer (no. 4 of tabled and 
an 0.86 Mq white dwarf for the point source and constant flux 
model (solid and dotted line), s' is compared to / for i] = 0.25 
and a = 0.1, 0.3, and 1.0 (short-dashed, dash-dotted, and dot- 
dash-dotted line, respectively). 



Figure 16. Mass transfer rate -M2 [Mq/yt] from an evolved 
initially 1.6 Mq MS star (no. 4 of tabled onto an initially 
0.6 Mq white dwarf for the point source model as function of 
orbital period P [hr]. System parameters: ij = 0.25, v = q, 
a = .3, and strong braking according to Verbunt & Z waanI 
(I1981ID . The dotted line shows the evolution without irradi- 
ation feedback. The onset of mass transfer cycles occurs at 
Ml * 0.91 Mq and M2 ~ 0.26 Mq. 



fig. shows^ that this is not the case, when the donor star 
reaches M2 = 0.45 Mq at an orbital period of ~ 8 hr. Instead, 
the system becomes unstable at a donor mass much lower than 
expected, i.e., at about 0.26 Mq. The reason for this discrep- 
ancy is probably that the bipolytrope model underestimates the 
timescale Tee on which the stellar radius changes if irradiation 
is changed. Numerical computations indicate that Tee from i66\ 
is too large by a factor of ~ 1.4 for a fully convective 0.3 Mq 
MS star whereas it is too smaU by ~ 1.6 for an 0.5 Mq MS 
star and even more so for more massive or more evolved stars. 
As a result evolved systems are significantly more stable than 
previously expected. 

6.1 .3. Giant donors 

Another group of CVs are systems containing a giant or sub- 
giant donor. Fig.^jshows the stability diagram for a CV with 
a giant donor of 0.8 Mq (model no. 6 of table Since Tee 
according to ( I66t gives an unreliable estimate for the thermal 
timescale of the convective envelope of giants, we have used 
the timescale Tce.u ~ 6 ■ 10^ yr which is defined as the ther- 
mal energy U of the convective envelope divided by the stellar 
luminosity. The timescale on which the giant radius changes, 
when irradiation is changed, is even larger than Tee.u by ~ 2. 
Giants typically have ^ ~ 10"^ and 6 ~ n + 5 ^ 2.ln our case 



* To perform the evolutionary computations shown in fies. ll6l and 
1221 from the initial turn-on of mass transfer to the minimum period 
without any discontinuity we used opacity tables for solar composi- 
tion in the outer layers. Comparisons with computations using correct 
Helium and CNO abundances for the final mass transfer phase indicate 
that, though using "incorrect" opacity tables in that evolutionary stage 
leads to a different stellar radius, it does not affect the mass transfer 
cycles significantly. 



for a - 0.3 the "island of instability" ranges from 0.25 to 20 in 
(/in ). This means the ratio of Tee and t^ must be in the range 



-2 



..5 ■ 10"^ < < 7 ■ 10 



or taking into account ( I89l l 



-2 



(94) 



(95) 



5.5 ■ 10"^ < — < 4 ■ 10 

Td 

Using analytical approximations given bv lRittejlll999t eq.43) 
and ( I42> we can write 



4M2 - Me / M2 \^ /0.25 M©^^ 



— ~-3-10 
Td M2 \MqI 



Me 



(96) 



for Pop. I giants and a similar expression for Pop. II giants. The 
temporal evolution of the core mass Me is given by 



Me = Me 



(97) 



Here, M^- i denotes the initial core mass and too the nuclear 
timescale, i.e., the time when the core mass formally becomes 
infinite. Almost all values on the right-hand side of ( I96> are less 
than unity. We have Me < M2 and for thermally stable systems 
typically M2 < Mq. Accordingly, is too small for destabi- 
lizing the system except if Me < 0.25 Mq and Me <sc M2. 
For the system shown in fig.[^Eq. J96> yields: 



5 ■ 10" 



Td 



(98) 



In fact, to destabilize the system a would have to be larger by 
about an order of magnitude than the value of 0.3 which we 
have chosen for this example. Since an efficiency a > 1 is 
rather unrealistic, we conclude that CVs with extended giants 
can not undergo mass transfer cycles. 
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Figure 17. s' for a CV with an 0.8 Mq giant (no. 6 of tabled 
and an 1 .2 Mq white dwarf for the point source and the con- 
stant flux model (solid and dotted line), s' is compared to / for 
7/ = and a = 0.1, 0.3, and 1.0 (short-dashed, dash-dotted, and 
dot-dash-dotted Une, respectively). 

There are three reasons, why CVs with giant donors are the 
more s table, the more extended the giants are: First, King et al. 
ill 9971) have neglected the last term in i79\ that restricts the re- 
gion of instability not only to high but also to low braking rates. 
Thus, systems do not get automatically unstable for sufficiently 
small braking rates. Second, the depth of the superadiabatic 
convection zone grows during giant evolution which shifts the 
"island of instability" to higher (/in). Hence, more evolved gi- 
ants are less affected by irradiation. Third, more evolved giants 
have a larger scale height ^ so that their "island of instability" 
is smaller than for less evolved giants. 

6.2. Low-mass X-ray binaries 
6.2.1. Unevolved MS donors 

As we have discussed in the last section, CVs might undergo 
mass transfer cycles for 0.1 < a < 1 given a suitable driving 
timscale. Another class of compact binary systems with ther- 
mally stable mass transfer, i.e., small q, are LMXBs. In the fol- 
lowing we consider only LMXBs with neutron star primaries. 
In principle, the discussion for these system is analogous to the 
discussion of CVs. While the radius of a white dwarf is about 
10"^ Rq, the radius of a neutron star is only about 10 km. Since 
for a given secondary star and driving mechanism (/in-) ~ 
(cf. Eq. ( I84» the smaller radius of a neutron star has to be com- 
pensated by a correspondingly small value of a for achieving 
a compara ble value of (/in). Because an efficienc y that small 
is unlikely, iKin g et alJ j 19961 Il997l) : iRitter et alJ GPOO) have 
concluded that irradiation in LMXBs is too strong to destabi- 
lize such systems. 

However, this conclusion appUes only to the constant flux 
model but not to the point source model. For the point source 
model s' decreases much more slowly for large (/in) than for 
the constant flux model because the surface elements near the 
terminator are only partially blocked by irradiation, even for 
very high fluxes due to the high angle of incidence. Therefore, 



the point source model does not require (/n) ~ 1 to get mass 
transfer cycles. Instead, s' becomes larger than / for fluxes be- 
low (/n) < 10^ - 10"^ as can be seen in fig.[T|lfor a LMXBs 
with an unevolved donor, in fig.^]for a LMXBs with a rem- 
nants of thermal timescale mass transfer, and in fig.|20|for a 
LMXBs with a giant. Thus, LMXBs can become unstable for 
a < 10-2 - 10-1. 
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Figure 18. s' for an LMXB with an unevolved 0.8 Mq MS star 
(no. 3 of tabled and an initially 1 .4 Mq NS for the point source 
and the constant flux model (solid and dotted line), s' is com- 
pared to / for ^ = 1 and a = 0.01, 0.1, and 1.0 (short-dashed, 
dash-dotted, and dot-dash-dotted Une, respectively). 
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Figure 19. s' for an LMXB with an 0.6 Mq remnant of ther- 
mal timescale mass transfer (no. 5 of tabled and an initially 
1 .4 Mq NS for the point source and the constant flux model 
(solid and dotted line), s' is compared to / for = 1 and 
a = 0.01, 0.1, and 1.0 (short-dashed, dash-dotted, and dot- 
dash-dotted line, respectively). 

The "island of instability" for the LMXB with an unevolved 
donor, which is shown in fig. ^] ranges from 0.02 to 150 in 
{fin-} for a = 0. 1 . This corresponds to 

6 ■ 10-^ < ^ < 4.5 ■ 10-^ (99) 
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Figure 20. s' for an LMXB with an 0.8 Mq giant (no. 6 of table 
[3 and an initially 1 .4 Mq NS for the point source and the con- 
stant flux model (solid and dotted line), s' is compared to / for 
77 = 1 and a - 0.01, 0.1, and 1.0 (short-dashed, dash-dotted, 
and dot-dash-dotted line, respectively). 



6.2.2. Highly evolved MS donors 

There is only a small region in the parameter space where 
LMXBs can undergo significant cycles for an irradiation effi- 
ciency as large as ~ 0.1 a nd for high braking rates according to 
IVerbunt & Zwaanl l ll98ll) . As an example we show in fig.ll9lthe 
evolution of a LMXB through thermal timescale mass trans- 
fer. The "island of instability" for the 0.6 Mq remnant (model 
no. 5 of table [3, corresponding to an orbital period of ~ 8 hr, 
is shown in fig.l22land ranges from 0.02 to 100 in (/„), corre- 
sponding to 

7 ■ 10"^ < ^ < 3.5 ■ 10^1 (101) 

With Tee ~ 5 ■ 10^ yr the driving timescale must be between 
2 ■ 10** yr and lO'^yr so that the system should be near the 
boundary of the unstable region. This is roughly consistent with 
the fact that the system becomes stable just one cycle later at 
P ~ 7 hr. 



With Tee ~ 2 ■ 10'' yr for an unevolved 0.8 M© MS star and 
using (I89> we get 

4.5 ■ 10** yr < Td < 3.5 ■ lO'^ yr. (100) 

From this we can infer that for ff = 0.1 such a system shows 
cycles for gravitational braking with Tgrav = y ~ 6 ■ 10^ yr and 
also for higher braking rates of up to 107giav. Numerical com- 
putations show that this system is rather near to the boundary 
of the unstable zone. Since the bipolytrope model underesti- 
mates Tee by a factor of ~ 2, the system undergoes cycles only 
for braking rates up to 57giav. Only low-amplitude sinusoidal 
oscillations occur for a while for a - 10"^ the system is 
deep in the unstable region as the evolution in fig. [^demon- 
strates. 
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Figure 21. Mass transfer rate -M2 [MQ/yr] from an evolved 
(Xc ~ 0.66) initially 0.8 Mq MS star onto an initially 1.4 M© 
neutron star for the point source model as function of orbital 
period P [hr]. Syste m parameters: ii - 1, a = 0.01, and strong 
braking according to Verbunt & Zwaan ( 1981). The dotted line 
shows the evolution without irradiation feedback. 
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Figure 22. Mass transfer rate -M2 [A^o/yr] from an evolved, 
initially 3.0 M© MS star (no. 5 of table Q onto an initially 
1.4M0 neutron star for the point source model and a = 0.1 
as function of orbital period P [hr]. For this computation fj - 
r] = I has been used as long as the Eddington accretion rate 
of the neutron star (2 ■ 10"^ Mo/yr) is not exceeded, otherwise 
fj - T] < 1. The dotted line shows the evolution of the same 
system without irradiation feedback. 



6.2.3. Giant donors 

Unlike CVs with giant donors LMXBs with giant donors can 
undergo mass transfer cycles even for very extended giants. 
This can be seen from the example shown in fig.|20|which illus- 
trates the conditions at the beginning of the evolution shown in 
Figs. |23] and [23 Unlike all other evolutions shown before this 
system evolves to longer orbital periods because it is driven by 
nuclear evolution, not loss of angular momentum. 

However, taking into account, that a neutron star has a quite 
small Eddington accretion rate of ~ 2 ■ 10"*^ MQ/yr, this par- 
ticular giant system would be stable, since the secular mass 
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Figure 23. Mass transfer rate -M2 [MQ/yr] from an 0.8 Mq 
giant (no. 6 of tabled onto an initially IAMq neutron star 
for the point source model, fj - I, and or = 0.1 as function of 
orbital period P [d] . The dotted Une shows the evolution of the 
same system without irradiation feedback. 
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Figure 24. Same as fig.|23lbut -M2 as function of time t [yr]. 

transfer rate is above the Eddington rate, so that there would be 
no irradiation feedback. 

6.3. General outline of mass transfer cycles and 
observational implications 

The appearance of mass transfer cycles differs between dif- 
ferent types of systems. "Ideal" outbursts as outlined by 
iKing etal.l ill 9971) require a very short ratio of Tee and which 
is typical for giants. They are characterised by a very short high 
state at a rather high mass transfer rate followed by a long last- 
ing low state as shown by fig.l24l 

In the limit of weak irradiation feedback, as it is typically 
the case for C Vs with 0.1 < a < 1 and weak braking, a system 
leaves the high state long before saturation is reached. Only a 
few 10 of the stellar mass is transferred during a cycle and 
the recurrence timescale is rather short. In the limit of strong 
irradiation feedback, as it is often the case for LMXBs with 
10"^ < a < 0.1, the system stays in the high state until satura- 
tion has been almost exactly reached. As a consequence of this. 



the mass transfer rate decreases to almost the stationary value 
before the system enters the low state. The duration of the high 
and low state are of the same order, and in the high state the 
system transfers mass just slightly above the stationary mass 
transfer rate for most of the time. 

Since the timescale of mass transfer cycles, i.e., the ther- 
mal timescale of the convective envelope of the donor is by far 
longer than typical observational timescales (at most of order 
of centuries), it is effectively not possible to observe the evo- 
lution of a single system from the high into the low state, or 
vice versa. Nevertheless, the occurrence of these cycles would 
affect the observable properties of a binary population. First, 
systems in the low state are more difficult to observe than in 
the bright high state or may even appear as detached systems 
if mass transfer ceases completely in the low state. If the dura- 
tion of the low state is longer than that of the high state, then 
most of the population can be basically unobservable. Second, 
even in the high state the mass transfer does not proceed at a 
constant rate but declines continuously from the peak rate until 
the system "switches" into the low state. Hence, even a com- 
pletely homogeneous sample can show significant variations in 
the mass transfer rate. Third, the mean mass transfer rate of the 
observable systems in the high state can be significantly larger 
than the secular mass transfer rate which can lead to apparent 
discrepancies between observations and theoretical predictions 
if the occurrence of mass transfer cycles is not taken into ac- 
count. 

Therefore, mass transfer cycles could provide an explana- 
tion for, e.g., the scatter of mass transfer rates of novalike CVs 
above the period gap, the disappearance of bright systems at 
the upper edge of the period gap, or discrepancies between pre- 
dicted and observed l uminosities and b irthrates of LMXBs as 
recently suggested bv lPfahl et al.l(l2003l) . 

6.4. Non-local effects 

We have assumed that irradiation is a local effect and can be 
treated for every surface element of the donor separately. If 
the contribution of non-local effects (e.g., circulations) to the 
lateral heat transport becomes non-negligible for sufficiently 
strong irradiation, or if there is another mechanism which heats 
the unilluminated side (e.g., scattering of X-rays by an ex- 
tended X-ray corona), then also the effective temperature and 
the intrinsic flux on the unilluminated surface depends on { fi„) 
and differs from Tq and Fo, respectively. 

Observations of the eclipsing binary AA Dor show that 
even in the case of very strong irradiation an illuminated star 
can preserve a cool hemisphere with ~ 2 000 K while the illu - 
minated side is heated up to ~ 20 000 K (iHilditch et alJl2003l) . 
Therefore, it seems plausible that also in the case of LMXBs 
at most a small fraction of the irradiating flux is transported 
to the unilluminated side. Nevertheless, for sufficiently strong 
irradiation, e.g., (/in-) ~ 10"*, even a small fraction can have 
a significant effect onto the unilluminated side. We note that 
with the general definition of s in ( I18> our analytical ansatz re- 
mains valid, even if the unilluminated side is heated up, and 
we formally get the same stability criterion. In this case the 
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shape of s' depends on the exact prescription of those non- 
local effects. Here we can not make any quantitative predic- 
tions. Nevertheless, it is plausible that s' becomes greater at 
higher fluxes ) if non-local effects are taken into account. 
This would mean that our local treatment of irradiation pro- 
vides a lower limit on the occurrence of mass transfer cycles in 
LMXBs. An upper limit is given by the fact that the area below 
the s' curve has to be less than the maximum fraction of the 
surface which is affected by irradiation, i.e., less than unity. 

7. Summary and Conclusion 

We have described the physics of irradiation driven mass trans- 
fer cycles by a two- dimensional system of ordinary differential 
equations similar to K ing et al ] (11996'). An additional term of 
order of several times ^ which reflects the usual thermal re- 
laxation of the mass losing star, can not be neglected, because 
it provides an upper limit for the driving timescale r^, respec- 
tively a lower limit for the timescale ratio where mass trans- 
fer cycles can occur Taking into account this term and also the 
rather large pressure scale height and the deep superadiabatic 
convection zone of giants we have concluded that CVs with gi- 
ants above a certain core mass can not undergo mass transfer 
cycles unless the efficiency parameter a is unexpectedly large. 
On the other hand LMXBs with giants might be suscept ible to 
the irra diation instability as it has been predicted by Kin ^ et alJ 
1I1997I) . 

We agree with results of iKing et alJ (Il996l) : iRitter et alJ 
regardin g CVs with unevolyed M S stars: For high brak- 
ing rates (e.g.. IVerbunt & Zwaanlll98ll), as they are also re- 
quired by t he period gap model (e.g.. ISoruit & Ritte3ll983t 
iKolbl 1 19931) . such systems are stable except for the most 
massive (~ 1 Mq) donor stars. On the other hand, for low 
braking rates, as they are proposed by Sills et alJ JioOO ): 
lAndronov et alJ (|2()03), CVs above the period gap might un- 
dergo cycles. In this case only novalikes are affected since 
dwarf novae have rather large amplitude outbursts with a rather 
small duty cycle and therefore a small a^. Even for gravita- 
tional braking the susceptability to irradiation ceases within the 
period gap so that CVs with unevolved MS stars below the pe- 
riod gap are stable. 

Numerical evolutionary computations have shown that the 
uncritical application of the predictions of the analytical model 
to highly evolved remnants of thermal timescale mass trans- 
fer does not yield good quantitative estimates for the bound- 
aries of the unstable region. The timescale on which the stellar 
radius changes, if irradiation is changed, seems to be signif- 
icantly underestimated by the bipolytrope model so that such 
systems are more stable at high braking rates than previously 
expected. The most evolved systems are susceptible to irradia- 
tion for 0.1 < a < 1 and orbital periods roughly between 2 and 
5 hours, or donor masses between 0. 1 and 0.3 Mq, respectively. 
For giants the ratio of thermal energy of the envelope divided 
by the stellar luminosity seems to give a good estimate for Tee. 

Since in the point source model surface elements near the 
terminator are irradiated at a high angle of incidence, they are 
only partially blocked by irradiation even for high fluxes. Thus, 
other than expected by .King et al.. a996. .1997.) : .Ritter et alJ 



teooob also LMXBs with unevolved or evolved MS stars might 
be susceptible to the irradiation instability. Regarding the driv- 
ing timescale for such systems basically the same restrictions 
apply as for CVs. LMXBs can undergo mass transfer cycles 
only for a < 0.1. The contribution of non-local effects of irra- 
diation, e.g., circulations, which were not taken into account by 
our model, might be non-negligible for LMXBs. It seems plau- 
sible that this effect tendentially destabilizes the mass transfer 
However, it is not clear in all cases whether mass transfer in 
LMXBs is driven by Roche lobe overflow or by irradiation- 
induced winds jBasko & Su nvaevI 1 19731 iBasko et alJ Il977t 
llben.Tr etalJll997l) ^ In the latter case our model could obvi- 
ously not be applied to these systems. 

We conclude that the analytical model gives a qualitatively 
and also quantitatively suitable description for the onset of 
mass transfer cycles except for highly evolved remnants of ther- 
mal timescale mass transfer We think that the boundaries of the 
unstable regions can be determined with an accuracy of a factor 
of < 2. 

To answer the final question whether the observable CV 
or LMXB population or parts of it is undergoing mass trans- 
fer cycles we need to know two basic parameters: the driving 
timescale, respectively the magnetic braking law, and the effi- 
ciency a. Nevertheless, we have shown that irradiation-driven 
mass transfer cycles in compact binaries are possible for not 
unreasonable values of a and t^. 
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